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EXPERIMENTAL  AND  THEORETICAL  DETERMINATION 
OF  Jic  FOR  2024-T351  ALUMINIUM  ALLOY 


by 


P.  W.  BEAVER,  M.  HELLER,  and  T.  V.  ROSE 


*■  SUMMARY 

The  i-integra!  is  an  elastictplastic  fracture  mechanics  parameter  which  can  be  regarded 
as  a  measure  of  the  intensity  of  the  crack  tip  stress  and  strain  fields,  irrespective  of  the 
plastic  zone  size.  The  value  of  )  at  the  onset  of  stable  crack  extension?)  has  been  sug¬ 
gested  as  a  fracture  criterion  for  both  large- and  small-scale  yielding  conditions.  In  this 
work  the  value  of  Jlc  for  an  extruded,  medium-strength  aluminium  alloy,  2024-T35!  bar. 
was  determined  using;  (i)  the  Hutchinson- Rice- Rosengren  crack  tip  model  and  experi¬ 
mentally-determined  crack  lip  strain  profiles ;  (ii)  a  finite  element-hybrid  contour  method, 
(iii)  a  modified  linear  elastic  fracture  mechanics  approach,  and  (iv)  the  ASTM  standard 
multiple-  and  single-specimen  techniques.  Agreement  between  the  values  obtained  from 
the  crack  tip  strain  profile  method,  the  two  numerical  methods  and  the  multiple-specimen 
methods  is  good. 

The  J  ,c  value  determined  by  the  single-specimen  method  is  not  valid  as  the  amount 
of  crack  extension  at  each  load  level  could  not  be  determined  to  the  accuracy  required  by 
the  ASTM  standard,  using  the  recommended  unloading  compliance  method. 
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1.  INTRODUCTION 


'A  The  application  of  linear  elastic  fracture  mechanics  (LEFM)  to  many  practical  situations 
is  inappropriate  because  crack  initiation  and  growth  are  usually  accompanied  by  crack  tip 
plasticity.  As  a  result,  considerable  effort  has  been  devoted  to  developing  methods  to  charac¬ 
terize  the  fracture  properties  of  metals  in  the  elastic-plastic  regime. 

The  .fefhtegral,  as  initially  proposed  by  Rice;[l  is  one  elasticyplastic  fracture  mechanics 
parameter  which  has  aroused  considerable  interest  in  recent  years.  It  can  be  regarded  as  a  measure 
of  the  intensity  of  the  crack  tip  plastic  stress  and  strain  fields,-similar  to  the  stress  intensity  factor. 
K,  in  LEFM.%egley  and  Landes  pJUrst  proposed  that  the  onset  of  crack  extension  under  plane 
strain  conditions  would  occur  when  the  .Aijitegra]  exceeds  a  critical  value  Since  then  consider¬ 
able  data  have  been  published  supporting  the  use  of  Jffas  an  elastic-plastic  fracture  criterion. 

This  paper  describes  a  test  program  to  determine  the  value  for  an  extruded  medium- 
strength  aluminium  alloy,  2024-T351,  which  is  used  in  the  aircraft  industry.  This  value  was 
determined  using:  (i)  the  Hutchinsonskice-Rosengren  crack  tip  model  and  experimentally- 
determined  crack  tip  strain  profiles,  (ii)  a  finite  element-hybrid  contour  method,  (iii)  a  modified 
linear  elastic  method,  and  (iv)  the  multiplec'and  single»tspecimen  methods,  as  described  in  the 
ASTM  standard  E-8 1 3-8 1 . 

V 


2.  BACKGROUND 

The  justification  for  using  ./-integral  as  a  ductile  fracture  criterion  is  based  on  the  Hut- 
chinson-Rice-Rosengren  crack  tip  model,  known  as  theTTRR  SlngulariTy  [3,  4).  These  authors 
determined  the  stress  and  strain  distributions  at  the  crack  tip  for  non-linear  elastic  materials 
which  display  a  Ramberg-Osgood  relationship  between  the  effective  stress  a  and  effective  plastic 
strain  c".  The  Ramberg-Osgood  relationship  is  given  by  the  following  equation: 

<T  =  <T,[cp]"  (I) 

where  n  is  the  work  hardening  exponent  and  <j,  is  a  constant.  McClintock  [5]  showed,  by  com¬ 
bining  the  HRR  model  with  Rice's  definition  of  J.  that  the  crack  tip  plastic  stress  and  strain 
fields  can  be  expressed  as  a  function  of  J  as  follows: 

<  =  <r,( 7/<r,  . (2) 

and 

el,  =  <-//<r,-/„),/ . . •£?,</»  (3) 

where  r  and  0  are  the  polar  co-ordinates,  /„  is  a  function  of  n  and  stress  state,  and  iu((l)  and 
elj (0)  are  universal  functions  of  0.  and  which  also  depend  on  n  and  stress  state.  The  exact  forms 
of  Equations  (2)  and  (3)  depend  on  the  stress-strain  law  used  to  represent  the  material  behaviour, 
and  whether  or  not  non-dimensionalizing  factors  are  used.  Consequently,  several  variations 
of  these  equations  have  been  published  in  the  literature,  for  example  see  References  6-8  In  this 
work,  the  forms  of  the  equations  as  derived  originally  by  McClintock  are  used. 

Comparison  of  the  plastic  stress  and  strain  equations  with  the  elastic  stress  field  equations 

from  LEFM  reveals  that  they  are  of  the  same  form.  Moreover,  for  the  case  when  n  =  I  (linear 

elastic  behaviour)  Equations  (2)  and  (3)  reduce  to  the  linear  elastic  fracture  mechanics  l/r1  3 
singularity. 


The  underlying  assumptions  used  in  deriving  the  HRR  equations,  however  place  some 
limitations  on  the  use  of  J  as  a  fracture  criterion.  Firstly,  the  material  behaviour  is  assumed  to 
conform  to  the  deformation  theory  of  plasticity,  that  is,  the  material  behaves  as  a  non-linear 
elastic  solid  and  no  unloading  occurs  during  deformation.  When  crack  extension  occurs  in  real 
materials,  some  unloading  takes  place  in  the  wake  of  the  new  crack  tip  [9],  Hence,  the  first 
limitation  is  that  J  can  only  be  used  to  represent  crack  initiation  under  monotonic  loading. 
However,  in  recent  years  J  has  been  used  for  predicting  fatigue  crack  growth  with  success  which 
suggests  that  this  limitation  may  be  relaxed,  for  example  see  References  10  and  II.  A  second 
consideration  in  using  J  is  that  the  region  ahead  of  the  crack  tip  in  which  J  describes  the  stress- 
strain  field  must  be  large  compared  with  the  microstructural  elements  of  deformation  and  frac¬ 
ture.  These  processes  occur  in  an  intensely-deformed  region  at  the  crack  tip,  the  size  of  which 
is  of  the  order  of  the  crack  tip  opening  displacement  (CTOD).  as  will  be  discussed  later.  There¬ 
fore,  the  radius  of  the  region  dominated  by  the  HRR  field  must  be  large  compared  with  the 
CTOD. 


3.  EXPERIMENTAL  PROCEDURE 

Aluminium  alloy  2024-T4  extruded  bars  of  rectangular  cross-section  (108  mm  x  |l  mm) 
were  rolled  to  6-4  mm  thickness  and  heat  treated  to  the  T35I  condition.  The  average  tensile 
properties  after  this  treatment  were:  0-2°o  proof  stress  426  MPa,  ultimate  stress  503  MPa 
and  total  elongation  17-2/„.  The  specimen  geometry  used  in  this  work  was  a  modified  compact 
tension  specimen,  as  shown  in  Fig.  I.  Numerically-determined  compliance  curves  have  shown 
that  this  specimen  is  more  compliant  at  shorter  crack  lengths  than  an  ASTM  compact  tension 
specimen  of  the  same  dimensions  [12]. 


3.1  ASTM  Standard  Multiple-  and  Single-Specimen  Test  Methods 

The  JIC  value  for  2024-T35I  was  determined  using  the  multiple-  and  single-specimen  test 
methods  as  specified  in  the  ASTM  standard  E8I3-8I.  Briefly,  these  methods  can  be  described 
as  follows: 

In  the  multiple-specimen  method,  several  specimens,  all  with  the  same  size  of  fatigue  crack, 
are  given  different  values  of  tensile  load  and  then  unloaded.  The  amount  of  stable  crack  extension. 
Aa,  for  each  specimen  is  determined  by  first  marking  the  crack  front,  and  then  completing  the 
fracture  statically  so  that  the  amount  of  crack  extension  can  be  observed.  (Crack  fronts  in 
aluminium  alloys  are  usually  marked  by  fatigue  cycling  at  low  loads  whereas  heat  tinting  is 
used  for  steels.  During  this  research  program  a  more  effective  method  of  marking  crack  fronts 
in  aluminium  alloys  was  developed.  This  new  method,  based  on  liquid  metal  embrittlement 
of  grain  boundaries  in  aluminium  by  gallium,  is  described  in  Reference  13).. 

The  single-specimen  technique  differs  from  the  multiple-specimen  technique  in  that  the 
Aa  values  are  estimated  from  the  unloading  slope  of  the  load  displacement  curve  and  the  elastic 
compliance  [14].  The  unloading  slope  is  determined  at  several  values  of  load  during  the  test. 
This  method  assumes  that  small  amounts  of  unloading,  i.e.  of  the  order  of  I0°„  of  the  load, 
will  not  disturb  the  fracture  process  but  will  provide  a  small  portion  of  a  linear  elastic  curve 
whose  slope  will  give  an  instantaneous  measure  of  the  crack  length.  The  compliance  calibration 
curve  used  in  this  work  is  given  in  Reference  12. 

The  values  of  J  for  each  Aa  value  in  both  techniques  arc  determined  using  the  area  under 
the  load  versus  load-line  displacement  curve,  A.  (Fig.  2(a))  and  the  following  equation:* 

J  =  Af(aa'W)!(Bb)  (4) 

where  f(aJW)  is  a  dimensionless  coefficient  which  is  a  function  of  the  initial  crack  length  a„ 
and  the  specimen  width  W.  B  is  the  specimen  thickness  and  b  is  the  remaining  ligament  width. 


*  Please  note  that  the  formulation  of  Equation  4  in  the  ASTM  standard  is  misleading, 
the  coefficient  f(aj h  )  should  be  in  the  numerator. 
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imen. 


DISPLACEMENT 
(ON  LOAD  LINE) 


a)  Several  specimens  are  loaded  to  different  values  and  the  corresponding 

J  values  are  determined  from  the  areas  under  the  load/displacement  curve. 


b)  A  J-resistance  curve  is  constructed  using  the  values  of  J  determined  in 
(a)  and  the  corresponding  crack  extension  (,'a)  values  determined  from 
the  fracture  surfaces. 


Fig.  2 


Schematic  diagram  of  the  procedure  for  determining  J ^  using 
the  ASTM  standard  multiple  specimen  test  method. 


A  ./-resistance  curve  is  then  constructed  by  plotting  the  values  of  /  for  each  load  level  as  a  func¬ 
tion  of  Aa,  as  shown  schematically  in  Fig.  2(b). 

A  crack  tip  blunting  line,  which  represents  an  amount  of  crack  extension  from  stretching 
of  the  crack  tip,  (as  distinct  from  fracture)  is  also  drawn  on  the  /-resistance  curve.  This  blunting 
line  is  approximated  by  the  following  equation 

/  =  2ay-  Aa  (5) 

where  ay  is  the  yield  strength.  In  addition.  0- 15  mm  and  1  5  mm  offsets  are  drawn  parallel  to 
the  blunting  line.  The  ASTM  standard  requires  that  at  least  4  data  points  fall  within  these  two 
offset  lines.  The  critical  JiV  value  can  then  be  determined  from  the  point  of  intersection  of  a 
linear  regression  line  fitted  to  these  data  and  the  hlunting  line.  This  point  represents  the  onset 
of  stable  crack  extension  from  the  blunt  crack  tip. 


3.2  Crack  Tip  Strain  Profile  Method 

J  can  also  be  determined  from  the  crack  tip  strain  distributions  at  various  points  on  the 
load  displacement  curves  using  a  photo-printed  square  grid  (40  lines  mm)  on  the  specimen 
surface  and  a  replicating  technique,  as  outlined  in  Reference  15.  Microgrid  patterns  were  applied 
to  the  surfaces  of  the  specimens  used  in  the  ASTM  single-specimen  tests.  The  variation  in  the 
spacing  of  grid  lines  was  less  than  0- 1  pm.  An  example  of  a  photo-printed  grid,  on  a  specimen 
which  has  been  fatigue  cracked,  is  shown  in  Fig.  3.  The  onset  of  stable  crack  extension  could 
also  be  determined  from  the  replicas,  and  hence  the  critical  load  for  first  stable  crack  extension 
/\  and  thc/K  value  can  be  determined,  as  will  be  described  later. 

The  values  of  /at  the  points  on  the  load  displacement  curve  at  which  replicas  were  taken 
were  determined  from  the  plastic  strain  profiles  ahead  of  the  crack  tip.  cpy{0  ~  0)  versus  r  data, 
using  a  75  |«m  gauge  length  and  Equation  (3).  An  example  of  a  typical  plastic  strain  profile  is 
shown  in  Fig.  4.  However,  strain  values  within  the  so-called  intensely  deformed  /one  [7]  were 
not  used  as  it  is  not  clear  whether  Equation  3  is  valid  in  this  region.  The  extent  of  the  intensely 
deformed  /one  in  front  of  the  crack  tip  is  approximately 


u  =  M.CTOD  (6) 

where  \f  is  a  constant  with  a  value  of  1  [ 1 6 j  or  2(7].  and  the  CTOD  is  given  by  the  following 
equation 


CTOD  -  NJ  ay  (7) 

where  N  is  a  constant.  Broek  (17)  found  that  .V  =  0-44  for  aluminium  alloys  so  that  at  the  critical 
load  CTOD  =  20  pm  which  is  within  the  experimentally  determined  range  of  values  using  the 
replicas.  In  the  present  work  no  strain  value  within  40pm  of  the  crack  lip  was  used  in  deter¬ 
mining  the  /  values.  Furthermore,  results  from  a  theoretical  analysis  of  Shih  and  German  [18] 
showed  that,  for  compact  tension  specimens,  the  HRR  singularity  is  only  valid  over  distances 
from  the  crack  tip  of  6  10  times  the  CTOD  value.  Therefore,  no  strain  value  greater  than  approxi¬ 
mately  200  pni  was  used  to  determine  the  /  values. 

The  values  of  /  at  each  load  level  were  determined  using  the  strain  values  within  the  above 
range  of  r  values  and  Equation  (3).  The  values  of  <t,  and  n  in  Equation  (3)  were  determined 
using  tensile  test  data  for  the  same  material  and  Equation  (I):  these  values  arc  678  MPa  and 
0  075  respectively.  Values  of  /„  and  f^(W)  are  given  in  Reference  '3;  for  a  tensile  crack  under 
plane  stress  conditions,  and  an  n  value  of  0  077  (the  nearest  tabulated  value  of  n  to  0  075). 
/„  =  2-87  and  -  0)  =  0-8.  (It  should  be  noted  that,  even  though  the  specimen  geometry 
satisfied  the  ASTM  requirements  for  plane  strain  conditions,  plane  stress  values  of  /n  and 
<^(fl  -  0)  were  appropriate  as  the  strain  profiles  arc  from  surface  measurements). 

The  values  of  /  for  the  same  load  displacement  points  at  which  the  replicas  were  taken 
were  also  determined  using  the  area  under  the  curve  as  for  the  multiple-specimen  technique. 
These  values  were  then  compared  with  the  values  determined  from  the  crack  tip  strain  profiles. 
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Fig. 


An  example  of  a  photo-printed  grid  used  to  measure  the  strains 
ahead  of  crack  tips.  This  photograph  is  of  the  fatigue  precrai 
J-integral  specimen  number  BP27  prior  to  loading.  The  lines  ai 
to  the  crack  are  intensive  slip  bands  produced  during  the  prec 


lastic  strain  profile  as  determined  froi 

rinted  square  grid  on  the  specimen  surfr 


Fig.  4  A  typical  crack  tip  p] 
replicas  of  a  photo-pi 


4.  NUMERICAL  ANALYSIS 


Two  numerical  methods,  both  assuming  plane  strain,  were  used  to  evaluate  J  for  the 
specimen  geometry  used  in  the  experimental  work.  Values  of  crack  length,  a  and  the  load  at 
initial  stable  crack  extension  P ,  were  obtained  from  the  experimental  results  for  the  single¬ 
specimen  tests;  for  specimen  number  BP22  a  =  25-55  mm  and  P ,  =  5-25  kN  and  for  specimen 
number  BP27  a  =  25- 10  mm  and  Pt  =  5  -  98  kN.  Both  numerical  approaches  used  the  finite 
element  method,  with  the  analyses  being  done  using  the  PAFEC  suite  of  programs.  The  stiffness 
matrices  were  computed  using  2x2  reduced  integration  and  double  precision,  and  the  solurion 
was  obtained  using  double  precision. 


4.1  Finite  Element — Hybrid  Contour  Method 

In  this  method  the  value  of  J  was  determined  using  the  displacement  results  of  an  elastic- 
plastic  finite  element  analysis  and  applying  the  hybrid  contour  method  [19].  The  finite  element 
mesh  used,  shown  in  Fig.  5(a),  consisted  of  212  eight-noded  iso-parametric  quadrilateral  elements 
and  11  six-noded  isoparametric  triangular  elements.  To  represent  the  loading  applied  during 
testing,  a  point  load,  P,  was  applied  at  the  top  of  the  pin  hole.  (It  has  been  demonstrated  in 
Reference  12  that  the  stresses  in  the  crack  tip  region  are  insensitive  to  whether  a  point  or  distri¬ 
buted  loading  is  applied  at  the  hole).  The  stress-strain  properties  of  the  material  were  approxi¬ 
mated  by  assuming  a  two-segment  piece-wise  linear  relation.  For  accurate  modelling  of  the 
plasticity  effects  surrounding  the  crack  tip,  special  crack  tip  elements  were  used,  and  the  total 
load  was  applied  in  increments.  The  relevant  details  are  as  follows: 


(i)  In  the  region  surrounding  the  crack  tip  the  mesh  is  particularly  refined,  with  the  two 
crack  tip  elements  being  approximately  I /60th  the  length  of  the  crack,  a.  These  two 
elements  had  their  midside  nodes  shifted  to  the  quarter  points  to  generate  the  near-tip 
r'12  displacement  singularity  [20],  Various  studies  have  shown  that  these  elements 
give  very  good  results  when  used  for  both  elastic  and  plastic  crack  tip  conditions  [21.  22). 

(ii)  The  first  load  increment  was  chosen  to,  approximately,  coincide  with  the  load  that 
first  caused  a  near  tip  gauss  point  to  become  plastic.  The  remaining  increments  each 
corresponded  to  5%  of  the  total  load.  A  solution  convergence  to  less  than  I  •3“„  error 
in  stresses  was  achieved. 


In  the  hybrid  contour  method  the  nodal  displacements  along  a  contour  in  the  elastic  material 
surrounding  the  crack  tip  plastic  zone  were  used  to  evaluate.  J  (These  nodal  displacements 
were  determined  in  the  finite  element  analysis).  The  contour  used  in  shown  in  Fig.  5(b)  along 
with  the  plastic  zone  as  also  determined  from  the  finite  element  analysis. 


4.2  Modified  Linear  Elastic  Fracture  Mechanics  Method 

In  this  approach  it  was  assumed  that  the  presence  of  crack  tip  plasticity  makes  the  specimen 
behave  as  if  it  contained  an  elastic  crack  larger  than  the  actual  size  of  the  crack,  a.  Hence  the 
effective  crack  size,  at„.  is  taken  to  equal  a  +  rFl 2,  where  rp  is  an  estimate  of  the  plastic  zone 
length  [23],  and  is  given  by 


K2 

■r  ( I  -68  a,)2 


(8) 


where  K  is  the  mode  /  stress  intensity  factor. 
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(a)  Finite  element  mesh  (half  of  specimen  modelled). 


Y 


(b)  Contour  used  to  compute  J. 


Fig.  5  Finite  element  modelling. 


An  iterative  procedure  was  used  to  estimate  the  value  of  K  at  the  crack  length  <r,„  for  both 
specimens  BP22  and  BP27.  Initially,  rr  was  determined  using  A  based  on  the  actual  crack  length , 
a  revised  K  was  then  estimated  for  a  crack  length  of  aef,  =  a  +  rP  2  The  various  values  of  A 
were  determined  from  the  displacements  (obtained  from  elastic  finite  element  analyses)  of  nodal 
points  behind  the  crack  tip  using  the  equation  given  in  Reference  24,  namely. 


uE  i2.Tr 
K  4(1  -  v2)  \J  r 


(9) 


where  u  is  the  load  direction  displacement  of  a  near  tip  node  on  the  crack  face  behind  the  crack 
tip,  r  is  the  distance  of  that  node  from  the  crack  tip  and  E  is  Young’s  modulus. 

Finally,  the  value  of  J  was  evaluated  using  the  estimated  value  of  A'  by  substituting  into 
the  equation  relating  J  and  K  for  elastic  plane  strain  conditions.*  namely. 

A2 

j  =  —(i  —  v2)  tio> 


where  v  is  Poisson's  ratio.  The  values  of  £  and  v  for  2024-T35I  are  £  —  72  4  GPa  and  i  =  0  31. 


5.  RESULTS  AND  DISCUSSION 
5.1  ASTM  Standard  Test  Methods 

A  Jlc  value  of  16-9  k/m2  was  first  determined  using  the  ./-resistance  curve  constructed 
from  the  multiple-specimen  test  data  as  shown  in  Fig.  6.  This  value  is  in  good  agreement  with 
the  estimated  Jlc  value  of  18-Ok.Am2  obtained  using  equation  10  and  the  A , ,  value  of  38  5 
MPas  m  for  this  alloy  as  given  in  Reference  25. 

The£K.  value  determined  from  the  single-specimen  ./-resistance  curve,  shown  in  Fig.  6.  is 
47-8  kV/m2.  This  value  is  much  higher  than  the  values  obtained  from  both  the  multiple-specimen 
test  data  and  the  estimated  •/,<  value.  The  large  difference  between  the  single-  and  multiple- 
specimen  J,c  values  is  unexpected  and  suggests  that  the  Aa  values  were  not  accurately  determined 
using  the  unloading  compliance.  This  was  verified  using  data  from  the  multiple-specimen  tests. 
The  Aa  values  for  specimens  used  in  these  tests  were  estimated  using  the  final  unloading  com¬ 
pliances  and  were  compared  with  the  values  measured  from  the  fracture  surfaces.  The  Aa  values 
determined  from  the  complaincc  values  were  at  least  27”„  less  than  those  measured  directly 
from  the  fracture  surface.  Therefore,  the  single-specimen  JH  value  is  not  valid  as  the  ASTM 
standard  requires  that  the  difference  between  the  estimated  and  real  Aa  values  must  be  less 
than  1 5° 

Three  factors  were  identified  as  being  responsible  for  the  large  difference  between  the 
estimated  and  real  Aa  values,  these  are: 

1.  non-linearities  in  the  mechanical  and  electrical  systems. 

2.  the  sensitivity  of  the  compliance  calibration,  and 

3.  crack  tunnelling  effects. 


*  It  should  be  noted  that  the  computer  program  used  to  determine  the  values  of  A  could 
simulate  clastic  plane  strain  conditions  even  though  the  specimens  were  only  6  4  mm  thick. 
(This  specimen  thickness  satisfies  the  ASTM  requirements  for  a  valid  plane  strain  test  whereas 
it  does  not  satisfy  the  ASTM  requirements  for  valid  linear  clastic  plane  strain  Au  test.) 
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Fig.  6  J-resistance  curve  for  2024-T351  as  determined  by  the  ASTM 
multiple  and  single  specimen  test  methods. 


The  first  factor  is  a  result  of  the  method  used  to  measure  the  unloading  compliance.  The 
recommended  method  is  to  electronically  subtract  the  elastic  part  of  the  load-displacement  curve 
and  replot  the  non-linear  part  with  a  higher  amplification.  This  amplification  improves  the 
accuracy  of  the  measurements;  however  small  non-linearities  in  the  mechanical  or  electrical 
system  become  apparent  during  the  unloadings.  For  example,  clip  gauge  hysteresis,  and  friction 
between  the  clip  gauge  and  knife  edges,  the  specimen  and  loading  pins,  and  the  loading  pins 
and  clevises  result  in  non-linearities  and  make  it  difficult  to  estimate  the  best  fit  lines  for  the 
unloading  curves.  The  errors  in  measuring  the  unloading  compliance  are  of  the  order  of  2",,. 
however,  the  resulting  errors  in  the  estimated  Aa  values  arc  up  to  100'.,,  i.e.  the  compliance 
calibration  is  very  sensitive  to  small  changes  in  the  unloading  compliance  value  The  sensitivity 
of  the  compliance  calibration  can  be  reduced  by  using  longer  initial  crack  lengths,  that  is  if  the 
a  W  ratio  is  increased  from  0*5  to  0-7.  where  If  is  the  specimen  width 

Examination  of  the  fracture  surfaces  alter  the  various  tests  showed  crack  tunnelling  which 
affects  the  unloading  compliance  and  hence  the  Aa  values.  When  the  crack  front  is  curved,  the 
crack  does  not  close  uniformly  during  unloading  as  it  is  not  closing  along  a  straight  hinge  line. 
This  uneven  closure  means  that  the  unloading  compliance  value  is  higher  as  the  shift  in  the  rota¬ 
tion  point  is  retarded  compared  with  the  situation  when  the  crack  front  is  straight  and  closure- 
uniform.  Therefore,  the  estimated  Aa  values  are  less  than  the  real  values  as  determined  by 
direct  measurement  (26)  and  this  results  in  the  7-integral  values  being  over-estimated.  Crack 
tunnelling  is  very  difficult  to  prevent  although  it  can  be  minimised  by  using  side-grooved  speci¬ 
mens  or  thicker  specimens. 


5.2  Crack  Tip  Strain  Profile  Method 

Crack  tip  strain  distributions  were  determined  from  replicas  of  the  surface  grid  patterns, 
examples  of  which  are  shown  in  Figs  7(a)  (d)  for  specimen  liP27.  The  plastic  strain  profiles 
ahead  of  the  crack  tip  vs  r)  determined  from  these  replicas  arc  shown  in  f  ig.  X.  I  his  figure 
shows  that  the  shape  of  the  plastic  strain  profile  changes  when  significant  crack  extension  occurs. 

The  values  of  7  at  each  load  level  were  determined  by  substituting  the  strain  values  at 
distances  of  between  approximately  40  and  200 /<m  from  the  crack  up  into  f  quatton  t5i.  These 
data  showed  that  the  several  values  of  7  at  each  load  level  were  essentially  the  same  1  he  mean 
values  of  7  at  each  load  level,  denoted  7STRMS.  and  the  standard  deviations  for  the  two  specimens 
used  for  this  technique  are  given  in  Table  I.  These  results  support  the  analytical  results  ol  Sh»h 
and  German  [18]  that  the  FI  HR  singularity  is  valid  for  distances  from  the  crack  tip  ol  up  to 
10  times  the  COD  value  for  compact  tension  specimens. 

The  values  of  7  at  each  load  level  were  also  determined  from  the  areas  under  the  load 
displacement  curves,  denoted  as  JXRl  A.  and  are  given  in  Tabic  I.  Agreement  between  thc7SIRAIN 
and  7ARhA  values  at  each  load  level  is  good  except  at  the  lower  load  levels  for  specimens  BP22. 

The  load  at  which  stable  crack  extension  first  occurred  and  the  associated  7(,  value  were 
established  by  examining  the  replicas  of  the  crack  tips  and  the  associated  grid  patterns  in  both 
an  optical  and  a  scanning  election  microscope.  As  mentioned  in  the  Introduction,  us  soon  as 
initial  crack  extension  occurs,  the  material  behind  the  new  crack  tip  starts  to  unload.  This 
phenomenon  can  he  seen  in  Figs  7(a)  (d).  As  the  load  is  increased  from  507  to  5  77  kN  the 
CTOD  and  the  width  of  the  secondary  crack,  region  A.  also  increased.  Figs  7(a)  and  (b).  How¬ 
ever,  as  the  load  is  further  increased  to  6  20  kN  the  CTOD  and  the  width  of  region  A  decreased 
significantly.  Figs  7(b)  and  (c).  This  change  in  the  crack  profile  is  a  result  of  the  unloading 
associated  with  initial  crack  extension  in  the  interior  of  the  specimen  (It  should  be  noted  that 
this  unloading  behind  the  crack  tip  did  not  relax  the  surface  strains  ahead  of  the  crack  tip.  as 
shown  by  the  strain  profiles,  in  Tig  8  corresponding  to  the  above  micrographs.)  A  small  amount 
of  crack  extension  also  occurred  at  the  surface.  Crack  extension  occurs  first  m  the  interior  of 
the  specimen  as  (he  triaxial  stress  state  increases  the  maximum  principal  stress  up  to  three  times 
the  uniaxial  yield  stress.  In  comparison,  the  maximum  principal  stress  at  the  surface  is  limited 
to  approximately  the  yield  stress.  On  further  increasing  the  load  from  6  20  to  6  50  kN.  Figs 
7(c)  and  (d).  significant  crack  extension  occurs  at  the  surface,  and  t!»c  CTOD  increases  although 
the  width  of  the  secondary  crack  decreases. 


II 


Fig.  7  Replicas  of  the  crack  tip  profiles  and  grid  distortions  for 
specimen  BP27  at  various  loads.  (a)  5.07  kN,  (b)  5.77  kN, 
(c)  6.20  kN,  and  (d)  6.70  kN  (elastic  limit  =  4.80  kN) . 


TABLE  1 


J-Iotegral  Values  at  Several  Points  on  the  Load  Versus  Load-Line  Displacement  Curves  as 
Determined  from  (a)  Crack  Tip  Profiles,  and  (b)  The  Area  Linder  the  Curve 


Specimen  Number 
and  Crack  Length 

Load 

(kN) 

/strain 

(ky/m2) 

Area  (,k«//rn2) 

Mean 

Standard 

Deviation 

Mean 

Standard 

Deviation 

BP22 

4-45 

10-5 

2-2 

14-4 

at  =  25  -55  mm 

4-88 

12-8 

2-5 

17-8 

0-7 

5-25 

20-0 

2-7 

21-4 

5-63 

=  70  ,um 

25-5 

1-0 

surface  crack 

growth 

BP27 

5-07 

19-7 

2-5 

18-0 

0-7 

a,  =  25  -10  mm 

5-77 

21  -2 

2-3 

22  1 

6-20 

22-4 

1-5 

25-0 

1  0 

S  40  (an 

34-3 

14 

surface  crack 

growth 

The  above  micrographs  show  that  initial  crack  extension  occurred  in  specimen  BP27  at  a 
load  between  5-77  kN  and  6-20  kN,  hence  Jlc  is  between  21  -2  and  22-4  kJ/m2.  An  estimate 
of  the  7,c  value  can  be  obtained  by  assuming  that  initial  crack  extension  occurred  at  a  load 
midway  between  the  above  values.  This  produces  a  J[c  value  for  specimen  BP27  of  about  21  -8 
kJ/m2.  Similarly,  the  replicas  showed  that  initial  crack  extension  in  specimen  BP22  occurred 
r*  a  load  between  5-25  kN  and  5-63  kN.  However,  the  J  value  at  the  higher  load  level  could 
not  be  determined  as  the  amount  of  crack  growth  at  this  load  level  significantly  altered  the  crack 
tip  strain  profile.  Therefore,  a  lower  bound  estimate  of  the  y,c  value  for  specimen  BP22  is 
20  0  ky/m2. 

The  accuracy  of  this  technique  obviously  depends  on  the  load  intervals  at  which  the  replicas 
are  taken.  The  best  procedures  to  estimate  J,c  are  to  either,  estimate  the  critical  load  using  the 
lower  bound  LEFM  estimate  of  y,c,  and  take  two  to  three  replicas  at  load  increments  approxi¬ 
mately  5%  on  either  side  of  P„  or  to  use  an  additional  specimen  to  experimentally  determine 
P,.  More  replicas  would  be  required  for  this  specimen  than  for  subsequent  specimens. 

As  mentioned  in  Section  3,  the  yK  values  determined  from  the  surface  crack  tip  strain 
profiles  are  for  plane  stress  conditions,  whereas  the  value  obtained  from  the  ASTM  multiple- 
test  method  is  for  plane  strain.  Under  elastic  conditions,  the  plane  strain  Jlc  value  is  less  than 
the  plane  stress  value  by  a  factor  of  (I  -  v2)  [16].  The  value  of  v  for  this  alloy  is  0-31  so  that 

y(c (plain  strain)  =  0-90,ylr  (plane  stress)  (II) 

This  equation  can  be  used  to  give  an  estimate  of  the  plane  strain  Jlc  values  using  the  plane  stress 
values  obtained  from  the  crack  tip  strain  profile  method  as  initial  crack  extension  was  associated 
with  only  a  small  amount  of  plastic  deformation.  These  plane  strain  Jlc  values  are  in  good 
agreement  with  the  J,c  value  obtained  from  the  multiple-specimen  test.  As  shown  in  Table  2(a) 
the  95%  confidence  intervals  for  the  two  methods  overlap  indicating  that  the  differences  between 
the  y,c  values  are  due  to  material  variability  and  experimental  scatter. 
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5.3  Numerical  Methods 


•Ac  values  were  also  determined  using  the  two  numerical  methods  outlined  in  Section  4. 
As  shown  in  Table  2(b),  the  Ac  values  determined  from  the  hybrid-contour  and  the  modified 
LEFM  methods  are  nearly  identical  for  each  specimen.  The  differences  between  the  values  for 
the  two  specimens  are  a  result  of  the  errors  associated  with  determining  the  load  P,  and  material 
variability. 


TABLE  2 

J)c  Values  Determined  By  Various  Methods 


Method 

Specimen 

Stress 

Ac 

95°, ,  Confidence 

Number 

State 

(kJim2 

Interval  (k/  mJ) 

(a)  Experimental 

Single  specimen  test 

plane  strain 

47-8 

(not  valid) 

Multiple  specimen  test 

— 

plane  strain 

16  9 

121-217 

Crack  tip  strain 

BP22 

plane  stress 

>200 

13-8-26  2 

Profile  method 

8P27 

plane  stress 

21  8 

16-2-27  4 

BP22 

plane  strain* 

2  180 

12  4-23  6 

BP27 

plane  strain 

19  6 

14-6-24-7 

(b)  Numerical 

FEM  hybrid-contour 

BP22 

plane  strain 

17  8 

method 

BP27 

plane  strain 

22-5 

14-5-26-1 

Modified  LEFM 

BP22 

plane  strain 

181 

approach 

BP27 

plane  strain 

22  5  J 

Using  KIC  from 

Ref.  25  &  cqn  10 

— 

plane  strain 

180 

*  Estimates  of  the  plane  strain  Ac  values  were  obtained  from  the  plane  stress  values  using 
equation  1 1 . 


The  determination  of  Ac  using  the  modified  LEFM  approach  is  based  on  the  assumption 
that  Equation  10,  which  comes  from  the  linear  elastic  definition  of/,  is  valid  in  the  clastic-plastic 
regime;  that  is,  if  the  amount  of  plastic  deformation  is  small  enough  Ac  will  be  identical  to 
G,c  and  hence  can  be  related  to  Klc,  However.  Ac  is  experimentally  determined  using  small 
specimens  which  reach  the  fracture  point  well  beyond  the  linear  elastic  regime.  In  comparison. 
Klc  is  experimentally  determined  using  thicker  specimens  which  reach  the  fracture  point  under 
linear  elastic  conditions.  This  means  that  the  measurement  point  for  Ac-  i.c.  at  initial  stable 
crack  extension,  may  not  be  coincident  with  the  measurement  point  for  A'u,  i.e.  at  2"„  crack 
extension,  especially  if  there  is  significant  plastic  deformation  [27],  Landes  and  Begley  [27)  used 
the  R-curve  to  illustrate  this  difference  in  measurement  point  between  the  two  methods.  For 
a  standard  ASTM  linear  elastic  AIC  test  the  R- curve  is  nearly  flat  so  that  the  A,,  value  at  2°„ 
extension  is  approximately  the  same  as  the  value  of  A'  at  initial  crack  growth,  as  shown  schemati¬ 
cally  in  Fig.  9.  Consequently,  Equation  10  is  valid  as  the  measurement  points  for  both  A,r  and 
Ac  are  comparable.  In  comparison  for  materials  which  fail  by  a  ductile  mechanism  or  for  thin 
sheet  “plane  stress"  behaviour,  the  R-curve  can  be  fairly  steep.  The  A',r  value  at  2”0  extension 
for  these  situations  would  not  be  the  same  as  the  K  taken  at  initial  crack  extension  Fig.  9.  This 
means  that  the  Ac  measured  at  the  point  of  initial  crack  growth  would  be  lower  than  the  A, 
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CRACK  AOVANCE.&. 


Fig.  9  Schematic  of  R-curve  showing  the  differences  uetween  the 
J1C  and  K1C  measurement  points  [27], 


value  determined  from  Aflc  at  2°0  crack  extension.  However,  in  this  work  the  load  at  which  first 
stable  crack  growth  occurred  P%,  as  determined  from  the  replicas,  was  used  to  calculate  A. 
denoted  as  K%,  and  not  the  load  at  2°„  extension  or  Penn-  Therefore,  the  values  of  J  determined 
from  the  Ky  values  using  the  modified  LEFM  approach  will  be  consistent  with  the  Ju  values 
determined  using  the  crack  tip  strain  profile  method.  Table  2  shows  that  the  Ju  values  obtained 
by  the  two  methods  are  in  agreement. 


6.  CONCLUSIONS 

The  Ju  value  for  extruded,  medium  strength  2024-T35I  aluminium  alloy  has  been  deter¬ 
mined  using  three  experimental  and  two  theoretical  methods  Agreement  between  the./,,  v;iluc' 
obtained  from  (i)  a  crack  tip  strain  profile  method.  ( 1 1 >  a  finite  element-h) hrid  contour  method. 
On)  a  modified  linear  elastic  farcture  mechanics  approach  and  (tv)  the  ASTM  multiple-specimen 
technique  is  good. 

The  JH  value  determined  by  the  single-specimen  technique  is  not  valid  as  the  amount  of 
crack  extension  at  each  load  level  could  not  be  determined  to  the  accuracy  required  by  the 
ASTM  standard,  using  the  recommended  unloading  compliance  method. 
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